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ABSTRACT 


A synthesis is made of several published studies of the normal heat budget of the lower troposphere over the 


Northern Hemisphere in winter. 


one using the thermodynamic energy equation and the other a heat-balance procedure. 


Charts of heat sources and sinks are presented, based on two independent methods, 


Although there are differ- 


ences in important details, both methods indicate that the horizontal scale of heating is the same as that of the normal 
lower-tropospheric temperature field, but that the field of heating is almost 90° out of phase with that of temperature. 


This result has an important bearing on the energy budget of the circulation. 


It is concluded that the heat- 


balance method is correct in indicating that there is a positive correlation between the normal temperature and 


heating fields. 


from heating to potential energy at the scale of the climatological long waves. 


maintaining these waves against friction. 


According to the energy transformation equations, this means that there is a positive transformation 


This energy is directly available for 


It is hoped that this study will be of some use in the design of numerical general circulation experiments, and 
as a basis of comparison with data from meteorological satellites. 


1. INTRODUCTION 


Attempts have been made by many investigators over 
the past 25 to 30 years to obtain an accurate picture of the 
distribution of sources and sinks of thermal energy which 
are responsible for driving the atmospheric circulation. 
Recent progress in improving the speed and capacity of 
electronic computers has raised the possibility of solving 
this problem, at least in principle, by including in a com- 
plex general-circulation model all the physical processes 
associated with heating and cooling. However, aside from 
the fact that these physical processes are poorly under- 
stood, there remains the difficulty that even the enormous 
speed and capacity of present machines are unequal to the 
task of resolving all of the complexities in a reasonable 
length of time. Therefore, in order to include some esti- 
Mate of the heating fields in general circulation experi- 
Ments, it has been necessary to employ simplified heating 


functions based largely on hypotheses arising from the 
judgment and experience of individual authors, as well as 
on the degree of complexity of their models (Phillips [25], 
Smagorinsky (to be published), Mintz [20], Charney [7], 
Fjgrtoft [9]). 

It is of vital importance that these estimates, though 
crude, give a reasonably accurate picture of the magnitude 
and pattern of the heat sources and sinks, because, as will 
be suggested at the conclusion of this article, their location 
relative to the different branches and space-time scales of 
the circulation is of critical importance in determining the 
sources of energy for, and therefore the life cycle of, these 
wind systems. Thus, there is evidently a need for a 
reasonably accurate picture of at least the normal fields of 
heating and cooling in the atmosphere so that these can be 
used as a check against the proposed heating functions. 

There is also some evidence that the fields of heating and 
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cooling are closely related to errors in current operational 
numerical circulation forecasts (Martin [19]). Therefore, 
a knowledge of how these fields are related to upper-level 
wave patterns should be of considerable practical assist- 
ance in routine forecasting. 

Finally, it is desirable to provide a rational basis of 
comparison with the directly observed values of the global 
heat budget obtained from meteorological satellites. 

For these three reasons the author was encouraged to 
attempt an up-to-date and, it is hoped, accurate synthesis 
of these older estimates of the normal heat sources and 
sinks in the lower troposphere over the Northern Hemi- 
sphere. Encouragement in this effort was provided 
through the recent publication of the results of two 
extensive studies on this subject (Staff Members, 
Academia Sinica [29], Budyko [5] [6}]). 


2. GENERAL PROCEDURE 


The first step was to make a search of the literature to 
locate all published studies of heating and cooling covering 
as much of the Northern Hemisphere as possible. This 
tusk was rendered fairly simple because of the relatively 
small number of projects treating the subject on so vast 
a geographical scale. Despite this restriction, it is 
possible that some of the pertinent literature may have 
been overlooked. Furthermore, the reader is cautioned 
that the list of references at the end of this paper is not 
meant to be an exhaustive bibliography on the subject, but 
contains only the reports actually used in this project. 

Next, the publications were divided into those in which 
the heating fields were determined using the thermody- 
namic energy equation, and those based on a_heat- 
balance method. It was hoped in this way to obtain two 
independent estimates, thus providing an accuracy check. 
However, it was soon discovered that the final results of 
the two methods were not comparable, and, as will be 
discussed later, it was necessary to make a choice as to 
which gives the more reliable results. 

After the pertinent reports for each of the two methods 
had been selected, the data were read from the published 
charts at each standard 10° intersection of latitude and 
longitude from the equator to 60° N., and converted to 
mean heating in the layer 1000 to 500 mb. expressed in 
c.g.s. units (eal. (It was found that 
not enough data were available north of 60° N. to obtain 
heating fields in the polar regions.) In only one case 
were original data tabulations used for extracting the 


gm.~! see.~! 10%). 


needed information. 

The effect of reading published charts, rather than 
copying data tabulations, is to decrease the magnitudes of 
the maxima and minima of heating. 

It was found that there were available no more than 
two independent sources of information dealing with any 
one particular type of heating field. These were simply 
averaged together in constructing the final charts, with 
no attempt to weight them in accordance with their 
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probable relative accuracy. To a certain extent such a 
screening for accuracy was taken care of in the original 
careful selection of source material. 

In some cases, only one source of data was available; 
and in one case, none at all. Thus, some reworking and 
extension of the original data was inevitable in order to 
obtain complete hemispheric heating fields. It must be 
stated that this was done in some instances by the use of 
approximate techniques which would probably not have 
been deemed entirely valid by the original authors, 
These modifications will be discussed in some detail below, 

The general forms of the equations used in the computa- 
tions are presented in Sections 3 and 4. These will not 
be derived here, as this has been taken care ot in some 
detail in the original source references. However, the 
major assumptions on which they are based, and _ the 
additional approximations used by the individual authors, 
will be treated here. 

The final heating fields presented in figures 1 to 5 may 
be considered to apply to the colder season and are 
averages for the layer 1,000 to 500 mb. or 0 to 5.5 km. 
No distinction was made between individual data fields 
for winter (December, January, and February), and those 
for the month of January only. Examination of normal 
flow patterns shows that there are relatively small changes 
in the principal meteorological parameters (and presum- 
ably also in the heating field) between the months of 
December and February. 


3. THERMODYNAMIC ENERGY EQUATION METHOD 


From the first law of thermodynamics and the definition 
of potential temperature, the individual rate of heating 
per unit mass, averaged over a specified time interval, 
may be written as follows: 
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Here, q is the amount of heat per unit mass; ¢,, specific 
heat at constant pressure; 7, absolute temperature; 1, 
time; V, horizontal wind vector; V, horizontal gradient; 
6, potential temperature; w, vertical wind component; 
and 2, height. The bar overscore represents the arith- 
metic-averaged value of a given quantity for the selected 
time interval, and the prime symbol denotes the departure 
of a quantity from its time-averaged value. 

The terms on the right of equation (1) may be desig- 
nated, in order, the local heating (heat storage term), 
horizontal heating (horizontal advection of heat), vertical 
heating (vertical advection of heat), horizontal eddy- 
advection of heat, and vertical eddy-advection of /eat. 

The two principal assumptions made in deriving equa 
tion (1) are relatively minor. They are: (a) The local 
rate-of-change and horizontal advection of pressure are 
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small compared to its vertical advection. (b) The quan- 
tity 7'/@ may be considered constant in time when used as 
a coefficient of the vertical heating term. 

The first attempt to compute average heating from 
equation (1) was made by Wexler [31]. The results re- 
poried on here are based on the more recent studies of 
Aubert and Winston [2] and Staff Members, Academia 
Sinica [29], who computed the normal heating in the 
lower part of the troposphere for individual months or 
seasons of the year. Both of these studies omitted the 
two eddy terms and used a constant or normal vertical 
stability in evaluating 00/02. They also used constant 
values of V, 7, and @ within a finite vertical layer of the 


atmosphere. This is equivalent to assuming that there is 


no correlation in the vertical between the horizontal wind 
components and temperature, or between vertical motion 


and stability. Neither source made computations north 
of about 70° N. Few estimates were made by Aubert and 
Winston south of 20° N., and by Staff Members south of 
10°. N. 

These two sources differ somewhat in other approxi- 
mations used in their work as follows: 

Aubert and Winston evaluated the heating for the layer 
sea level to 10,000 feet using normal sea level and 10,000- 
ft. charts. The gradient wind obtained from the average 
pressure field for this layer was used in computing V. 
The quantity 7 is the average mean virtual temperature 
between these two levels, and ® was obtained from the 
divergence of the mean gradient winds. Their results 
must be questionable near high mountains, because pres- 
sures reduced to sea level were used, and no mountain- 
effect was included in the computation of @. 

The Staff Members used normal constant-pressure charts 
to evaluate heating in the layer 1,000 to 500 mb. (approxi- 
mately sea level to 5.5 km.). The geostrophic wind at 
700 mb. was used in computing V, and 7 was replaced 
by the thickness between the two layers. The mean 
vertical motion was obtained using the average divergence 
at the ground and at 500 mb. The authors state that the 
surface divergence was computed from observed winds, 
although it is difficult to see how this is possible in most 
areas of the globe, where no winds are available. The 
upper-level divergence was computed with the aid of the 
vertically integrated tendency equation. This seems to 
avoid the difficulty arising from the balance of terms in 
the divergence equation. The effect of topography was 
included in computing vertical motion at the ground, but 
this apparent improvement of the method of Aubert and 
Winston is largely nullified by the use of reduced 1000-mb. 
charts in computing thickness advection. Thus, their 
heating fields are probably also unreliable in mountain 
areas 

In order to obtain an average of these two heating fields, 
it Was necessary to convert the values of mean heating in 
the liver 0 to 3 km., derived from data of Aubert and 
Winston, to the layer 0 to 5.5 km. This was done using 
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the rather arbitrary assumption that the intensity of 
heating per unit mass varies linearly with height, reaching 
a maximum at 1 km. and vanishing at sea level and 10 km. 
Using this scheme, the ratio between average heating in 
the two layers can easily be obtained. Thus, the 
conversion was accomplished simply by dividing the values 
obtained from Aubert and Winston by 1.08. 

The two independent estimates are remarkably similar 
(separate charts not shown). The major difference is in 
eastern Asia and the western Pacific, where the heating 
maximum computed by Staff Members is weaker and 
considerably farther to the east than that of Aubert and 
Winston. This discrepancy is difficult to trace, but may be 
due to the use of different normal circulation patterns. 

The average of these two heating fields is shown in 
figure 1. A brief discussion of the pattern of heating and 
its relation to the broadscale temperature field will be 
deferred until Section 5. Here, an attempt will be made 
to assess its probable accuracy. This is not an 
matter, in view of the numerous assumptions discussed 
previously. 

Inspection of figure 1 quickly reveals certain regions 
where one would suspect that the values are in error. 
For reasons given below, the indicated heating over eastern 


easy 


Canada and Siberia appears to be too large in comparison 
with that in the northwestern Atlantic and Pacific Oceans. 
Also, there is no apparent explanation for the small but 
positive heating over the oceans west of California and 
North Africa, which incidentally, appears on the charts of 
both authors. The northeastern continents and south- 
eastern oceans in winter are known to be regions where 
precipitation is relatively light. This means that con- 
densation, the only important source of heat in these 
areas, must be relatively weak. On the other hand, in both 
regions the air sometimes loses heat by contact with a 
colder land or ocean surface, and always loses heat by 
radiation to space. 

Examination of the individual terms in equation (1) 
reveals the probable source of these errors. The reliability 
of the normal circulation patterns is probably quite 
acceptable, so that, except in the high mountain areas, 
the sum of the first two terms is reasonably correct. Thus 
the errors may be traced to the omission of the eddy terms 
and to the approximations used in computing vertical 
motion. 

The necessity of accurate estimates of vertical motion 
stems from the fact that there is a strong negative corre- 
lation between the sum of the first two terms and the 
third term on the right side of equation (1). This has 
been demonstrated on a synoptic scale by Panofsky [24] 
who showed that vertical motions computed by the 
adiabatic method (essentially equation (1) with zero heat- 
ing and the eddy terms omitted) are positively correlated 
with weather and with vertical motions computed from 
the divergence of observed winds (kinematic method). 
Several recent extensive studies of the kinetic energy 
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Figure 1.—-Normal winter heating, sea level to 5.5 km., 
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budget of the circulation have been based on vertical mo- 
tions computed by the adiabatic method, evidently with 
some assurance that this gives reasonably accurate results, 
at least for that purpose (e.g., Jensen [15]). 

Aubert and Winston [2] have shown that for nearly 
time-averaged circulations there is also a 
between the horizontal and 


steady-state 

high negative correlation 

vertical heating terms. 
This tendency for a negative correlation between three 


computed by the thermodynamic energy equation method. 
*day~! multiply by 4; to C °. day~!, by 4/100. 
+ 40 units cross-hatched; cooling greater than - 
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of the important terms means that equation (1) must be 
used with extreme caution in estimating heating, since a 
large systematic error in any term may even lead to the 
wrong sign of the heating. This may result in serious 
errors in the spatial pattern as well as the magnitude of 
heating. This precaution applies especially to the term 
containing vertical motion, because all the known (indi- 
rect) methods of estimating the vertical wind component 
are subject to systematic errors. 
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in the two studies discussed here there is evidence that 
the computed vertical motions are too small. Indeed, 
Aubert and Winston show that for a few selected localities 
the mean monthly divergence (and therefore vertical mo- 
tion) computed from observed winds is four times as large 
as the divergence of the gradient wind. This means that 
in general there is not enough compensation between the 
horizontal and the vertical terms. In regions of gen- 
erally sinking motion in the northeastern continents and 
southeastern oceans this results in too much heating, while 
in the areas of slow mean upward motion in the western 
oceans, not enough heating is computed. 

The omission of the eddy terms also has quite serious 
consequences in certain areas. Attempts to estimate the 
horizontal eddy effect (fourth term on the right of equa- 
tion (1)) by an approximate theory of horizontal mixing 
have been made by Moller [21] and Elliott and Smith [8]. 
If it can be assumed that horizontal mixing is propor- 
tional to the horizontal temperature gradient, with a 
constant coefficient of proportionality, then it can be 
shown that: 


Vv =—AvT 2) 


where A is called the ““Austausch” coefficient (always posi- 
tive according to theory), and ¥?7 is the horizontal 
Laplacian of the mean temperature field (often called the 
thermal vorticity). 

In regions of strong cyclonic thermal vorticity, as in 
the northeastern continental areas (see fig. 6), the eddy 
term is negative and, like the vertical motion term, tends 
to reduce the magnitude of heating given by the horizontal 
term alone. An opposite effect is produced in regions of 
anticyclonic thermal vorticity, as indicated by anti- 
cyclonic shear of the isotherms in the western oceans 
(fig. 6) and anticyclonic curvature in the northeastern 
oceans, ‘. 

An attempt was made in the present study to include 
quantitative estimates of the eddy term as computed by 
equation (2), but this was abandoned in view of the 
obvious crudeness of this method as revealed by the 
studies of Méller and Elliott and Smith. 

It is clear that the only way to avoid this problem 
completely is to compute the terms using synoptic maps 
and then average them over the desired time interval. 
This task, too laborious to accomplish on a large scale by 
the older “hand” methods, has recently been tackled 
with the aid of an electronic computer by Wiin-Nielsen 
and Brown (to be published). Their preliminary results 
for one winter month show that a region of weak cooling 
replaces the heating in figure 1 in eastern Canada. How- 
ever, this may be due to the fact that the circulation for 
the chosen month (January 1959) was extremely abnormal 
at hich latitudes. 

The conclusions to be reached from the above arguments 
are that the strong heating centers in the northeastern 
continents in figure 1 should be weakened and shifted to 


the southeast, and that the weak heating in the south- 
eastern oceans should be eliminated, or even replaced by 
centers of cooling. 


4. HEAT-BALANCE METHOD 


The sum of all sources of real heat gain or loss in a 
certain column of air of unit cross-section is equal to the 
flux-divergence plus the storage of heat in that column. 
With the aid of a continuity equation for real heat, it can 
be shown that the storage plus flux-divergence is equal to 
the total individual rate of heating. Thus, the average 
heating per unit mass in the column is given by the 
equation: 

f-£ (r-+R+P). (3) 
Here g is the acceleration of gravity; Ap, the pressure 
difference between bottom and top of the column (here 
500 mb.); Zr, the net heat gain due to condensation or 
evaporation; ?, heating due to radiation; and P, that 
due to turbulent exchange with the earth’s surface. As 
before, the bar overscore represents an average over a 
specified interval of time. 

The only approximation in this equation is to omit 
other sources of heating or cooling which are generally 
considered negligible. Among these sources is the trans- 
formation between kinetic energy and heat in the fric- 
tional boundary layer. 

As a matter of practical convenience, it is assumed in 
addition that (a) the net condensation and evaporation 
is measured by the heat equivalent of rain or snow reach- 
ing the ground; and (b) the radiation balance is deter- 
mined by the net radiative inflow at the bottom and top 
of the column. In accordance with assumption (a) the 
quantity Z in equation (3) is the latent heat of condensa- 
tion or sublimation, and r is the amount of precipitation 
reaching the ground. 

In practice the individual physical processes on the 
right of equation (3) are evaluated from climatological 
data. Thus, depending on the equations used, a certain 
residual effect is present, which is due to correlation in 
time between parameters appearing in these equations. 
The consequences of these “eddy” terms will be treated 
briefly in a later paragraph. 
that this is probably not very important. 

In addition to being relatively free of errors associated 
with eddy terms, the heat-balance method has another 
important advantage over the thermodynamic energy 
equation method in that there is no tendency for com- 
pensation between individual terms, in the sense that 


Here it will only be stated 


there is a high negative pattern correlation between any 
two of them. In fact, just the opposite is true, as will be 
clarified in later paragraphs. 

To compensate for these advantages, it is difficult to 
formulate the individual physical processes in mathe- 


matical terms. Even if this were possible, the necessary 
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Figure 2.—Normal winter heating, sea level to 5.5 km. due to conc 
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hemispheric data for evaluating them are frequently 
lacking. 


A. EVAPORATION AND CONDENSATION 


To evaluate this component of heating according to 
previous assumptions, it is necessary to determine the 
hemispheric distribution of winter precipitation, and 
simply convert this to heating units by multiplying it by 
the latent heat of condensation or sublimation. The 
result is then substituted for Lr in equation (3) under the 
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assumption that all condensation takes place below the 
500-mb. level. 

Two separate fields were assembled, one using precipi- 
tation over the land by Landsberg [16] and over the oceans 
by Jacobs [14], and the second using hemispheric pre- 
cipitation patterns worked up by Moller [22]. These were 
averaged and multiplied by a constant latent heat of 
condensation, equal to 600 cal. gm.~! The results are 
shown in figure 2. Strictly speaking the latent heat of 
sublimation should have been used in those areas were 
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precipitation reached the ground in the form of snow. In 
addition, the magnitude of latent heat should be a func- 
tion of the temperature of condensation or sublimation. 
However, these refinements would make little difference 
in the magnitude and pattern of heating. 

Another minor discrepancy arises from the assumption 
that the net condensation and evaporation is measured by 
precipitation reaching the ground. This is a good assump- 
tion except possibly in those areas where clouds form and 
drift into other regions before evaporating or discharging 
their moisture. However, the heat content of such clouds 
is probably quite small, and it is difficult to isolate large 
areas Where a systematic effect of this kind is to be found. 

More important than these considerations is the ques- 
tion of the accuracy of the precipitation amounts. Some 
encouragement is gained by noting that the two precipi- 
tation fields are very similar, both in pattern and magni- 
tude (separate charts not shown). However, this must 
be tempered by the realization that the authors used 
similar procedures and sources of data. 

Over the lands, precipitation amount is determined with 
fair reliability from direct observation. 

Over the oceans, where direct observations are lacking, 
precipitation amount is determined from the frequency of 
rainfall. This is converted to precipitation by relating 
it to known amounts at island and coastal stations, after 
taking into account the effect of latitude and geographical 
location. The net result is to cast doubt on the magni- 
tude, but not the pattern, of rainfall over ocean areas. 

Two other sources of error influence the computed 
magnitude of heating, especially in the western oceans. 
The first is the truncation effect, mentioned earlier, due to 
reading precipitation charts at widely spaced grid points, 
This results in a reduction in the amplitude, or difference 
between maxima and minima. This is partly compen- 
sated by the second error, resulting from the fact that in 
regions of heavy rainfall cloud tops frequently extend 
above the 500-mb. level. Thus, part of the condensation 
occurs above the 1000- to 500-mb. layer, resulting in an 
overestimate of heating in these regions. 

If the above conclusions are accurate, namely, that at 
least the pattern of heating shown in figure 2 is correct, 
this helps to establish the pattern of the total heat budget 
of the troposphere. For, as pointed out by Albrecht [1], 
and as will be shown further on in this report, there is a 
high positive correlation between the precipitation pattern 
and that of total heating. 


B. SENSIBLE HEAT EXCHANGE 


The gain or loss of real (sensible) heat by direct tur- 
bulent contact with the earth’s surface was obtained for 
ocean areas from estimates by Jacobs [14] and Budyko 
[5]. Jacobs made use of the so-called “Bowen formula” 
which is based on an assumed simple ratio between the 
turbulent exchange of heat and that of water vapor. 
The evaporation is assumed to be given by the product 
of the wind at observation level times the vertical vapor- 
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pressure gradient. The sensible heat exchange is thus 
a function of evaporation, temperature gradient, and 
water vapor gradient near the sea surface. 

Budyko used a formula giving the heat exchange as the 
product of the wind and the vertical temperature gradient 
multiplied by a certain exchange coefficient which is 
assumed to be a constant. 

These formulae are regarded by many authorities as 
being fairly reliable, in spite of the undeveloped state of 
boundary-layer turbulence theory. The exchange co- 
efficient, assumed constant, is known to depend to some 
extent on stability. Normal winter values of the different 


parameters, determined with reasonable accuracy from 
ships’ observations, were used in evaluating the equa- 


tions. 

The two separate heating patterns (not shown) are in 
good agreement, although it must be stated that essen- 
tially the same sources of oceanic data were used. 

No winter charts of sensible heat exchange over land 
ure known to the author. However, Budyko [5] presents 
a chart for the year as a whole, as well as values for each 
month of the year for 12 stations representing certain 
selected climatic zones. These values are obtained from 
the equation for heat balance of the land surface, express- 
ing the sensible heat exchange as the difference between 
radiation of all wavelengths reaching the ground and the 
heat of evaporation. The amount of evaporation is in 
turn determined as the difference between total precipita- 
tion and runoff. 

In order to obtain a rough estimate of the pattern of 
sensible heat exchange over land areas, the present author 
adopted the following scheme: (a) Using Képpen’s cli- 
matic classifications, the land areas were subdivided into 
climatic zones approximating the definitions used by 
Budyko. (b) Within each zone, the difference between 
January heating and that for the year as a whole was as- 
sumed to be the same as that for the appropriate climatic 
station for that zone. (c) Budyko’s yearly maps were 
then read at standard latitude-longitude intersections, 
and the results of step (b) applied to these yearly means. 
In order to read Budyko’s maps, the lines of constant sen- 
sible heat exchange were extrapolated through mountain- 
ous areas although Budyko was careful to discontinue the 
lines in the higher elevations. 

The results of this admittedly crude procedure are 
plotted in figure 3 along with the average of the two ocean 
Only the general magnitude and sign of the 
This shows, 


estimates. 
computations over land have any reliability. 
as one might expect, that the values are in general negative 
over northern continents (where relatively warm air flows 
over a cold snow surface) and positive over the southern 
continents (where cold air from the north flows over 
warmer snow-free ground), and that the largest magnitude 
(—16 units) is considerably smaller than those over the 
oceans. The general magnitude of heat loss over the 
northern continents is in fair agreement with independent 


estimates of others. For example, Hanson [12] has re- 
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FIGuRE 3. 


cently computed a value of —11 units at the South Pole 
in spring. 

Comparison of figures 2 and 3 shows a high positive 
correlation between these two independent heating fields, 
as mentioned previously. This correlation is especially 
good north of 20° N. 

C. THE RADIATION BALANCE 


The third important component of the heat-balance is 
obtained by subtracting the net loss of radiation by all 


wavelengths at the bottom of the atmosphere from the 








Normal winter heating, sea level to 5.5 km. due to exchange of sensible heat with the earth’s surface. 
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See figure 1 for legend. 


net incoming radiation at the top. This radiation flux- 
divergence for the entire atmosphere is then “reduced” 
to the layer 1000-500 mb. 

The net incoming radiation is obtained from Simpson's 
[28] classical computations, which appear to remain the 
only available source for the geographical distribution of 
this quantity. According to Godske et al. [11], Simpson's 
charts probably portray the pattern of radiation quite 
well, although recent advances in computation met! ods 
might well improve the magnitudes. 
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simpson first computes the effective incoming short- 
wave solar radiation, and subtracts from this the net out- 
going long-wave radiation. In his calculations, both of 
these quantities depend largely on cloudiness, latitude, 
and season. The outgoing long-wave radiation is obtained 
by Simpson’s well-known technique of first dividing the 
black-body radiation curve into finite regions depending 
on the transparent or opaque parts of the water vapor 
spectrum. The outgoing radiation to space in the trans- 
parent regions depends on the normal surface temperature, 
while that in the opaque regions is assumed to depend on 
the temperature of the uppermost water vapor layer, 
located in the lower stratosphere. The temperature of 
the stratosphere is taken as a function of latitude only. 
When clouds are present, they are assumed to radiate as 
a black-body at constant temperature (260° A.). Simp- 
son’s method appears to result in too small values for the 
outgoing radiation, since later investigators have found 
that the effective upper water vapor layer lies well below 
the tropopause. 

The net radiation loss of the air near the ground is taken 
from Budyko [5]*. This depends on a rather complicated 
chain of four empirical formulae of the type pioneered by 
Angstrém and Brunt. The formulae depend on clima- 
tological mean values of albedo, cloud amount, the fourth 
power of surface air temperature, vapor pressure in the 
lower atmosphere, and the fraction of short-wave radiation 
reaching the ground with cloudy skies. The albedo 
depends on the latitude, season, and average type of 
ground cover; while the fraction of short-wave radiation 
reaching the ground with cloudy skies is assumed to be a 
function of latitude only. 

In extracting data from Budyko’s radiation maps the 
isolines were extended through mountain areas as in the 
case of the sensible heat. Also, after the results had been 
subtracted from Simpson’s values and then analyzed, 
certain discontinuities depicted by Budyko along coast- 
lines and the shores of inland water bodies were ignored. 

This net radiative flux-divergence for the whole atmos- 
phere was reduced to the layer 1000 to 500 mb. with the 
aid of London’s [17] mean zonal cross-sections showing 
the local radiative heat loss as a function of altitude and 
latitude for each season. From London’s cross-section 
for winter, a coefficient was obtained (for each latitude) 
equal to the ratio between the mean heating in the layer 
1000-500 mb. and that for the entire layer 1000-0 mb. 
This was multiplied by the total flux-divergence for the 
appropriate latitude to obtain a reduced value. Finally, 
the latter was adjusted so that its zonal mean for each 
latitude agrees with those of London. 

The final result of the various steps is shown in figure 4. 
Unlike the field of figure 3 the pattern correlates very 
poorly (although still slightly positively) with the net 


°A ilar study has been made by Bernhardt and Phillips [3]. Unfortunately, their 
comp! te radiative heat balance, including the long-wave components, was not available. 
Howe. or, it may be noted bere that their charts of total short-wave radiation reaching 
the ¢ 1 agree very well with those of Budyko. 
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precipitation pattern shown in figure 2. It is easy to show 
qualitatively that this poor correlation is due largely to 
the method of reduction. This may be demonstrated 
with the aid of two schematic diagrams prepared by 
Moller [23] showing the variation with altitude of radiative 
flux-divergence with clear skies, and with cloudiness and 
rain. It is clear from these diagrams that in regions of 
small precipitation, where there are frequent clear skies 
or low clouds having tops well below the 500-mb. level, 
over half of the total radiative flux-divergence lies below 
500 mb. Therefore in these regions (central and northern 
continents) the magnitude of radiative loss in the lower 
troposphere should be somewhat larger than that computed 
with a reduction factor depending on latitude alone. On 
the other hand, in regions of heavy precipitation, (western 
and northern oceans and near the equator) where cloud 
tops frequently reach to or above the 500-mb. level, the 
radiative flux-divergence shown in figure 4 must be 
sharply reduced because it is almost zero within and below 
the clouds. These modifications would clearly lead to an 
increase in the positive correlation between the fields of 
radiative heating and condensation. 

No attempt has been made in this study to apply such 
modifications in a quantitative fashion. It was felt that 
this would be tantamount to “begging the question.’ 
The solution to this problem must await a fresh approach 
to the computation of radiative flux-divergence as a 
function of geography and elevation. 


, 


D. TOTAL HEATING 


The sum of the three principal components of heat- 
balance is shown in figure 5. This is the final estimate 
of the normal heating of the lower troposphere by the 
heat-balance method. As expected, it shows a high 
positive pattern-correlation with the distribution of 
heating by precipitation (fig. 2) despite the unsatisfactory 
results for the radiative component (fig. 4). The lack of 
any major modification of the pattern is due to the fact 
that the radiation distribution is almost zonal, with 
relatively weak gradients. Correction of the field of 
radiation would have the effect of increasing the magnitude 
of the maxima and minima in figure 5 (especially the 
maxima) without changing the pattern greatly. 

Albrecht [1] has computed the field of heating for the 
year as a whole, also using the heat-balance method. His 
pattern is quite similar to that in figure 5 over the oceans 
and along the eastern coasts of continents. As might 
be expected from seasonal changes in heating over land, 
his chart shows less cooling in continental interiors. 


E. THE EDDY TERMS 


It must be clear from the discussion in subsections A 
through D above that the most important question regard- 
ing the heat-balance method is that of the validity of the 
crude semi-empirical procedures used for evaluating the 
components. For this reason it seems almost beside the 
point to consider the possibility of additional errors due 
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Figure 4. 


Normal winter heating, sea level to 5.5 km. due to radiation of all wavelengths. 
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See figure 1 for legend. Relative maxima 


underlined. 


to correlation between the individual meteorological! 
variables used in evaluating the equations. Furthermore 
the formulae for radiation and sensible heating have been 
designed specifically for use with climatic data, and would 
probably not be used if these two heating components 
were to be evaluated on a day-to-day basis. Neverthe- 
less, the problem of eddy effects is so important for the 
thermodynamic energy equation method that it will also 
be raised in connection with the heat-balance method. 


Therefore, a brief discussion of its possible importance 
will be given here. 

There is no eddy effect in the computation of heating 
by evaporation and condensation because this is obtained 
simply by summing the daily amounts or frequencies of 
precipitation. Therefore, it may be concluded that the 
pattern of precipitation (and presumably also the paitern 
of total heating) is free of errors from this effect. 

The formula used by Budyko for sensible heat exchange 
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lay be subject to some error due to a positive correlation 
between the exchange coefficient and the temperature 
gradient. This arises because the turbulence exchange 
coeflicient becomes smaller the more stable is the lapse 
rate. This would lead to some increase in the estimate of 
heating over ocean areas portrayed in figure 3. However, 
Budvko [6] states that this effect is not important over 
oceans. 


The crude estimates of this sensible heat component 





Normal winter heating, sea level to 5.5 km., computed by the heat-balance method. 
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See figure 1 for legend. 


over land would fortunately not be further biased by eddy 
effects because the heating is estimated by linear heat- 
balance equations. 

Since Jacobs used the Bowen ratio, the effect of a 
variable coefficient is partly eliminated. The good agree- 
ment between Jacobs’ and Budyko’s values over oceans 
provides some weak evidence that eddy effects are 
unimportant. 

In the radiation computations there are several oppor- 
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Fiaure 6. 


Normal January thickness, 1,000 to 500 mb., after Jacobs [13]. 


May 1961 


Lines drawn for every 100 meters and labeled in meters. 


Centers of minimum thickness labeled C; maximum, W. 


tunities for correlations to arise between parameters such 
as the albedo and short-wave radiation reaching the 
ground. Budyko [6] has attempted to make estimates 
of the errors arising in various steps of the radiation 
budget at the ground. He concludes that there is a 5 to 
10 percent error in the monthly values of short-wave 
radiation, but that larger errors arise in the total budget 
because of irregular variations in albedo and = surface 
temperature. However, in comparing estimates of total 
radiation with direct instrumental observations for a few 


Russian stations, he gets what he terms “very satis- 
factory” agreement. His computations also agree with 
the independent calculations by others. 

Thus we may conclude that eddy effects play a rela- 
tively minor role in the heat-balance method. 


5. THE ENERGY SOURCE FOR THE CIRCULATION 


The two independent heating patterns (figs. 1 and 5) are 
quite dissimilar in their details, but their large-scale fea- 
tures are qualitatively alike, and show a similar overall 
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Figure 7.—Average normal winter heating and thickness in the 
laver 1,000 to 500 mb. and in the latitude band 30° N. to 60° N. 
as a function of longitude. Solid curve, heating by the thermo- 
dynamie energy equation method. Dashed curve, heating by 
the heat-balance method. Dash-dot curve, thickness departure 
from latitude average. 


relationship to the mean temperature field in the lower 
troposphere, portrayed here by the normal thickness be- 
tween 1000 and 500 mb. (fig. 6). The latter has been con- 
structed from normal data recently published by I. Jacobs 
[13]. 

Both heating patterns show a broad ring of maximum 
heating near the equator, associated with the highest 
tropospheric temperatures; maxima in middle latitudes in 
and to the east of the thermal troughs along the coasts of 
Asia and North America, and in the central Mediter- 
ranean; and cooling in and to the east of the long-wave 
ridges in western North America and the eastern Atlantic. 
Both charts show a broad region of cooling over the 
Eurasian continent associated with a rather flat west-to- 
east thermal pattern. 

These large-scale relationships are brought out more 
clearly at middle latitudes in figure 7, which shows the 
average thickness and heating in the belt 30°N.to 60°N. 
plotted as a function of longitude. The curve of heating 
by the thermodynamic energy equation method agrees in 
almost all details with an independent computation by 
Wiin-Nielsen [32]. 

Figure 7 also emphasizes the considerable difference in 
magnitudes of heating and cooling by the two methods 
and shows an important phase difference between them. 
The maxima and minima of heating by the thermodynamic 
method are displaced 10° to 20° of longitude farther to the 
west. This seemingly small phase difference is quite im- 


portant when considering the energy budget of the circu- 


lation. Lorenz [18] has shown that in order to get a posi- 
live transformation from thermal to potential energy, 


Which is available for driving the circulation, there must 
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Figure 8.—Zonal averages of normal winter heating and thickness 
in the layer 1,000 to 500 mb. as a function of the sine of latitude. 
Solid curve, heating by the thermodynamic energy equation 
method. Dashed curve, heating by the heat-balance method. 
Dash-dot curve, thickness. Latitudes of maxima and minima of 
heating indicated. 


be a positive correlation between heating and temperature. 
But the correlation between these two quantities, as por- 
trayed in figure 7, is negative (—0.42) according to the 
thermodynamic and positive (+ 0.39) according to the 
heat-balance method. Thus, depending on which method 
is assumed correct, one may conclude that the semiperma- 
nent circulation in middle latitudes in winter is either 
energy-consuming or energy-producing. From the pre- 
vious arguments in Sections 3 and 4 it is probable that the 
maximum heating centers of the thermodynamic method 
are too far west, while the pattern obtained from the heat- 
balance method is essentially correct. Thus, it is con- 
cluded that the semipermanent systems are energy-pro- 
ducing on the average. Some further interesting details 
of the energy budget will be given in later paragraphs. 

Figure 8 shows the zonal averages of heating and thick- 
ness plotted against a sine function of latitude in order 
that a unit distance on the abscissa may represent a fixed 
area on the earth. Both of the heating curves show a 
positive correlation between thickness and heating south 
of 20° N. and north of 45° N. According to the equations 
of Lorenz [18], this means that in these regions there is a 
positive transformation between zonal heating and zonal 
available potential energy. The cooling at 20° N. and 
heating at the equator, shown in figure 8, may be identified 
with the sink and source of heat for the Hadley Cell, the 
thermally directed meridional circulation of low latitudes. 
The results of this study indicate that there is a positive 
transformation between heating and zonal potential 
energy at these latitudes, providing an important source 
of energy for the Hadley Cell. This is in agreement with 
a study by Tucker [30] which indicates that within the 
Hadley Cell there is also a positive transformation between 
zonal potential and zonal kinetic energy. 
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Figure 9.—-Normal winter heating as a function of the sine of 
latitude after Pisharoty [27], solid curve; and after Gabites [10], 
dashed curve. See text for explanation. 


It is possible to provide a rough check on the validity 
of the latitudinal heating curves by comparing them with 
other studies. Two of these are summarized in figure 9, 
although unfortunately neither one has been computed 
for the same atmospheric layer, period of time, or heating 
components as those in figure 8. The reader must be 
referred to the sources for a discussion of the assumptions 
and approximations of these two estimates. 

The first estimate, shown as the solid curve in figure 9, 
represents the horizontal geostrophiec flux-divergence of 
heat in the layer 1000-200 mb. and for January and 
February 1949 (Pisharoty [26]). This omits the flux- 
divergence of heat by the meridional cells. The second 
estimate, shown as the dashed curve in figure 9, gives the 
normal heat-balance for Jaruary for the entire earth- 
atmosphere system (Gabites [10]). In addition to the 
average individual heating for the entire atmosphere, this 
includes heating in the oceans, measured by the latitudinal 
flux-divergence of heat produced by ocean currents. 

Comparing these curves with those in figure 8, we note 
that there is good overall agreement in the location of the 
maxima and minima of heating. The curves of figure 9 
are in better agreement with the magnitude of the heating 
maximum in middle latitudes obtained by the thermo- 
dynamic method, and suggest that the corresponding 
“maximum” for the heat-balance method is much too low. 
On the other hand, they are in better agreement with the 
location (phase) of the mid-latitude heating “maximum” 
for the heat-balance method. Both of these observations 
are in agreement with conclusions reached earlier from an 
analysis of the two procedures. 

Finally, it is of some interest to attempt a more detailed 
quantitative analysis of the thermal energy budget of the 
semipermanent systems using the heating fields computed 
by the heat-balance method. The reader is referred to 
the original paper of Lorenz [18] for derivation of the 
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Figure 10.—Production of eddy available potential energy by the 
covariance of normal mean temperature and normal heating for 
winter, plotted as a function of the sine of latitude. 


functions. Here, a somewhat 
different derivation by A. Wiin-Nielsen (unpublished 
lecture notes) has been used. This shows that the trans- 
formation between total heating and total available 
potential energy per unit area is given by the equation: 


energy transformation 


Ap 
~ STq 


J 7+ drdy. (4) 


Here, C is a measure of the production of total available 
potential energy by the semipermanent systems (energy 
from this source which is available for maintaining the 
total mean circulation against friction); Ap, the pressure 
thickness of a thin vertical layer of the atmosphere; T and 
“4, the absolute temperature and individual heating in 
the layer; and S, the area of the earth’s surface for which 
the mean energy transformation is to be computed. The 
asterisk represents departure of a quantity from its area 
average, and the bar overscore, as before, represents an 
average in time. 

Following Lorenz, one can divide the total energy pro- 
duction into two components: One is the production of 
eddy available potential energy (energy available for 
maintaining the semipermanent long waves); and the 
other the production of zonal available potential energy 
(energy available for maintaining the north-south tem- 
perature gradient). The production of eddy available 


; , , — dq* 
energy is obtained simply by replacing 7* and = by the 


departure of each quantity from its zonal average; while 
each of these two quantities is replaced by the departure 
of its zonal average from its area average in order to com- 
pute the production of zonal available energy. 
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TasLeE 1.—Computation of total, eddy, and zonal available potential 
energy production per unit area in two selected latitude belts of the 
Vorthern Hemisphere for winter. (Units are ergs cm.~? sec.—) 
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In this study it has been assumed that the mean tem- 
perature and heating do not vary in the vertical through 
the 1000- to 500-mb. layer, so that Ap is equal to 500 mb. 
Also, the mean temperature in the denominator of the 
coefficient has been set equal to a constant. 

Figure 10 was computed from equation (4) using the 
data of figures 5 and 6, and shows the production of eddy 
available potential energy per unit area as a function of 
latitude. This suggests that at low latitudes the semi- 
permanent long waves provide only a small fraction of the 
energy necessary for their own maintenance. Most of 
this energy must come from other sources, such as that 
generated by eddies of smaller time scale. However, 
north of 40° N., the energy produced by the semi- 
permanent waves increases rapidly, until at 60° N. it is 
the same order of magnitude as that removed by friction. 
The latter is either 2000 or 4000 ergs cm.~* sec.~', de- 
pending on whether one accepts the estimates of Pisharoty 
(26] or of Brunt [4]. 

The total energy production (per unit area) and its two 
components have been computed for the entire area of the 
Northern Hemisphere between the equator and 60° N., 
and in the mid-latitude westerly belt between 30° N. and 
60° N. The results are summarized in table 1. 

Qualitatively, this table shows results which are to 
be expected from other considerations: in the region of 
the Hadley Cell and other meridional cells, a relatively 
large amount of zonal energy is produced. This is the 
bulk of the total energy production for the whole area 
and especially at low latitudes, while the relative im- 
portance of the zonal energy production diminishes 
greatly at higher latitudes. Perhaps somewhat surprising 
is the result that more total energy is produced per unit 
area at high than at low latitudes. 

The relatively small magnitudes in table 1 are of some 
concern. The largest is only about 10 to 20 percent of 
the average frictional dissipation. This suggests that 
even at high latitudes the semipermanent circulations 
must be supported by energy from systems of smaller 
time scale, perhaps that produced in the traveling cyclones 
and anticyclones. However, these figures are probably 
much too small for two reasons: the magnitudes of the 
heating maxima at middle latitudes are too small, as 
concluded in this and previous sections; and there is a 
ver) important positive correlation between heating and 
temperature in the vertical, contrary to one of the 
assumptions made in evaluating equation (4). The 
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vertical distribution of temperature and heating is 
important for the semipermanent long waves as well as 
for the meridional cells. This must be especially true 
for the Hadley Cell, where all the heat of condensation is 
released in the lower troposphere in the Doldrum belt, 
while the principal heat sink is at the top of a strong 
moisture discontinuity in the upper troposphere. If 
these two corrections could be included, they would 
result in a marked increase in the covariance between 
heating and temperature when integrated both hori- 
zontally and vertically over the selected region 


6. CONCLUSIONS 


There are differences in important details of the two 
normal heating patterns (figs. 1 and 5), which indicates 
that either or both of them are in error. However, they 
agree that the phase difference between the normal 
heating and temperature fields is close to 90°. There- 
fore, any procedure used in computing the heat budget 
must at least yield a good estimate of the pattern of heat- 
ing, because a small error in pattern can change the sign 
of the energy transformation. Also, it must be clear that 
theoretical heating functions used in proposed general 
circulation models must also adhere to this critical phase 
relationship if they are to be at all realistic. 

Future estimates of the thermal budget should be made 
by a fresh attack using the heat-balance approach, and 
they should include the vertical as well as the horizontal 
distribution of heating. Sources and sinks of energy for 
the summer season and for the Southern Hemisphere 
should also be considered. 
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ABSTRACT 


The development and motion of fronts associated with lee-side troughs on large mountain barriers has been 
investigated. These fronts differ from ordinary cold fronts in their horizontal temperature field, which is charac- 
terized by a sinusoidal thermal ridge. The thermal ridge intensifies, while remaining stationary with respect to the 
mountains, and moves eastward upon the approach of a Pacific cold front. 

An equation is derived, showing that changes in the thermal pattern can be described by changes in a potential 
thermal vorticity equation, which consists of three terms: (1) one representing an advection of the potential thermal 
vorticity by the 500-mb. wind; (2) one representing the advection of 500-mb. absolute vorticity by the thermal 
wind; and (3) a purely orographic term. 

An idealized sinusoidal model of the thickness pattern is used in conjunction with the prognostic equation to 
explain the development and motion of lee-side thermal ridges. Actual examples from synoptic maps are chosen 
to corroborate the theory. The conclusions are: (1) the thermal ridge will develop when the surface flow is such 
as to produce large-scale descent on the lee slopes of the mountains; (2) no thermal ridging will appear when the 
500-mb. ridge lies east of the lee slopes; (3) thermal ridging will appear with the approach of a 500-mb. ridge from 
the west; and (4) the thermal ridge will move eastward upon the passage of the 500-mb. ridge. 


1. INTRODUCTION 


A well-known phenomenon is the lee-side trough, asso- 
ciated with large mountain barriers. At certain times of 
the year low-pressure troughs frequently develop on the 
eastern slopes of the Rocky Mountains. A particular 
type of front can be found along the axis of the trough, 
which is distinguishable from ordinary cold or warm 
fronts by its horizontal temperature field. This kind of ~/fotherms 
front will be called a “‘pseudo front” in this paper. Here 
the word “pseudo” is not to be confused with other defi- 
nitions involving the term, as applied to squall lines and 
other meso-features [1]. Typically, cold fronts lie along 
the axes of low-pressure troughs, moving from colder to 
warmer air. In a typical cold front the thermal gradient 
is quite pronounced behind the front (fig. 1). Pseudo 
fronts differ from cold fronts in that (1) the associated 
surface trough does not show an abrupt temperature dis- 
continuity at low levels, although the wind shift across 
the front may be quite pronounced, and (2) the thermal 
contours are almost symmetrical about the surface trough 
axis, With a maximum temperature at the front, decreas- 
ing at about the same rate on either side of the trough. 
It is evident that this description is exactly that of the 
temperature field associated with an occluded front. 
Pseudo fronts are distinguished from occlusions on the 





Cold Front 


Pseudo Front 











‘Based on a Master’s Thesis in the Department of Meteorology, Massachusetts 
Institute of Technology, and part of a project prepared for the Air Force Cambridge s 
Rese ch Laboratories, Contract AF (604)-5491. Figure 1. 


Schematic diagrams of cold and pseudo fronts. 
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Figure 2.—-Cross-section, Winnemucca, Nev. (WMC), to Topeka, 

Kans. (TOP), 1200 Gur, January 13, 1959. 

therms of potential temperature; dashed lines, vertical velocity, 
(12 hr.)-'. Shading indicates greater than three-eighths 
Distance is in miles. 


Solid lines are iso- 


mb. 
cloud cover. 


basis of the development of the latter in the life of a 
frontal wave. In the early stages of cyclogenesis, the 
thermal structure of a pseudo front often resembles that 
of an occlusion, except that the system is intensifying. 
There seems to be no evidence in the historical continuity 
of pseudo fronts that a cold front has overtaken a warm 
front, such as occurs toward the end of the evolutionary 
evele of a classical Bergeron cyclone. 

The most common type of pseudo front on the North 
American continent develops along the lee slopes of the 
Rocky Mountains, particularly during the months of 
November through March. The front first appears in the 
thermal pattern as a sinusoidal ridge, with a wavelength 
of about 1500 km. At the surface, a marked wind shift 
line appears coincident with the thermal ridge axis. The 
wave intensifies in both the thermal and surface circula- 
tions, often becoming extended along the entire eastern 
slopes of the mountains, although during initial develop- 
ment the forward advance is often very slight. McClain 
{2} and Hess and Wagner [3] have studied in detail the 
development of lee-side cyclogenesis, and their analyses 
of 12-hour synoptic maps showed that the lee cyclones 
began to move eastward upon the approach of a trough 
from the Pacific. Cloudiness and precipitation accom- 
panying the front are noticeably absent [4]. A cross- 
section through both the approaching cold front from the 
Pacific and the pseudo front at 1200 amt, January 13, 1959, 
reveals not only the total lack of middle cloudiness and 
precipitation accompanying the pseudo front, but points 
out its distinct identity in the temperature field (fig. 2). 
The phenomenon of the lee-side pressure trough is of such 
frequent occurrence in the Rocky Mountain region, that 
mean sea level charts show it clearly [3]. The time re- 
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quired for full development of a lee-side trough is about 
36 hours, and the rapid intensification of the thermal ridge 
east of the Rocky Mountains is associated with the com- 
monly experienced “chinook” or “foehn.”” The ultimate 
eastward movement of the pseudo front seems to be influ- 
enced by approaching Pacific troughs. 

The purpose of this paper is to try to explain the devel- 
opment and movement of these pseudo fronts in terms of 
the large-scale physical processes in the atmosphere. 
Similar work has been done by McClain [2], in which he 
correlated upper tropospheric divergence with the develop- 
ment of the lee-side troughs. Other researchers such as 
Colson [5] and Newton [6] have tried to explain the fre- 
quency of cyclogenesis in the central United States in 
terms of air flow over the Rocky Mountains. 


2. THEORY 


Since pseudo fronts are best revealed in the low-level 
temperature pattern, it becomes necessary to discuss the 
motion of the accompanying thermal ridge. It would be 
best to derive a comprehensive prediction equation for 
the thickness pattern, from the basic laws governing the 
atmosphere. One approach, which works reasonably well 
and is simple, is use of a two-parameter model which 
describes the atmosphere in terms of the mean wind and 
a thermal wind, taken over an appropriate depth in the 
atmosphere. The following modeling assumption is 
chosen: 

Viz, y, P, O=Vanls, Ys OF BDV ola, y, OD (1) 


where V refers to the non-divergent component of the 
wind, and B is a function of pressure only. The definition 
of the mean quantities is written 


: 1 "Po 
( | = — | 
Po Jo 


The subscripts m and 0 respectively refer to the mean and 


( )dp. (2a) 


1000-mb. levels, where 


( )mo=( )m—( )o- (2b) 
From the previous definitions, it is evident that B,,—0 
and that the vertical shear is constant in a particular air 
column. Applying the curl operator to (1) gives the 

corollary 
5=Smt Bono, (3) 


where (is the relative vorticity component in the vertical. 
Introduction of the simplified vorticity equation 
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where V is the gradient operator at constant pressure, 
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w= A and f is the Coriolis parameter, yields an expres- 
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sion for large-scale atmospheric motion. The simplifica- is st 9 
: : ee ; a= ’ (12) 
tion refers to the deletion of the twisting, vertical advec- Jap 


tion, and frictional terms, and the assumption that f>¢. 
The boundary conditions for the model atmosphere are: 


w(p=0)=0, (5a) 
w(p= ps) =o, (5b) 


where subscript s refers to the ground surface. Sub- 
stitution of (1) and (3) in (4) yields 


of Om _ 


P+ BSF = — Vin WS tf) — BV im + Vm 


— BV yn 0 En +f) — BM « Vim toe (6) 


” ° ° ra) m a 
lo eliminate oe, apply (2) to (4), obtaining another ex- 


. ra) ry’ 
pression for So and subtract the result from (6). Thus, 


ine 7 —Vino* Vint D—{ Vat Vno[ = =]}. *Vomo 


f t“~ a 
¥5 B (= Bp (7) 
Application of this equation at the level where B’= (B*),, 


and defining B, = y(B*),, yields: 


OFmo 
ot 


aii f Ws 
+ Be 
(8) 


==Vno “Vn tf) Ven Vina ($2) 


Terms containing w can be eliminated by introduction of 
the thermodynamic considerations. For a temperature 
field where the potential temperature 6=6(z,y,p,t), 


dé 086 00 
a . oA 9 
di att ¥ V4 © Op (9) 
The quantity @ is defined in the usual manner where 
=T (™) —% Pi 
o=T (7) R po? (10) 


where Tis the absolute temperature, a the specific volume, 
R the gas constant for dry air, pp>=1000 mb., and «=R/e,. 
" O00\ 8 (Oa . 

Since (3) ~ a) and V@=6/aVa, equation (9) 
becomes, after ceabaatieatinns by a/@ and substitution of 
these relationships 


ad _ Oa 
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The requirement that the atmosphere be in hydrostatic 
equilibrium gives 


where g is the gravitational constant and z is the height 
of a constant pressure surface. Substitution of (12) into 
(11) yields, upon division by g and rearrangement of terms 


st (s)-— i v(5)+ 


Assuming that no non-adiabatic heating occurs, integra- 
tion of (13) with respect to pressure from 1000 mb. to the 
mean level, where the pressure is p», yields 


a O00 a dé 


+50 op” godt (13) 


ma OO 


O80 F » 
x ot — Vv, Vemat ("SX dp, (14) 
, Oz Oz Zz . 
since V,,5°V (3,)=9 and V- v(5,)=V= -V (35) From 
the equation of continuity 
Pp Pp *p 
w=—f V+ (Vat BV mo) dp=! wT -Veo | Bap. 
(16) 
Similarly 
Ow Ws - 
(35), =< S— Ba “no (17) 


dD 
By substitution of (16) in (14), and defining c= | Bdp the 
0 


thermodynamic equation becomes 
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ry Vin Vemot 7 H Ve Vno D, (18) 
where Em [ree -dp, and D= 3 x dp. Substitu- 


tion of (17) into (8) yields 
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Elimination of ¥-Vno between (18) and (19), and re- 
arrangement of terms yields 
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; 
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which can be written in the abbreviated form 


Vino ' VitntS)—kew,, 
(21) 


a a — =_— . —kK- am 
ot (¢ m) K2Zmo) — Vin V( Smo K2 no) 


where K=f/D and k=fE/(p.D). E and D can be evalu- 
ated from Buch’s [7] mean data, and from mean values of 
06/Op for western North America. 
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In a typical synoptic situation, changes in {mo and 
—Ke2no augment each other since increases or decreases 
in thermal vorticity are associated, respectively, with 
decreases or increases in thickness. The combined expres- 
sion ({mo—AZno) can be called a potential thermal vor- 
ticity. This means that the thermal vorticity no has a 
potential value with respect to some value Of 2m. Using 
the symbol @ for the potential thermal vorticity, expres- 
sion (21) becomes 


3) 


4 -—V,, . Vbo—V no ~ VitntS) —kw,, 


(22) 
dg 


i (23) 


—Vm0»V(Sm+f)— hey. 

The terms —Vao- Vien+f) and —ke, represent non- 
conservative, or generative, effects. The presence of 
generative terms suggests a means by which thermal ridges 
are suddenly built up on the lee slopes of large mountain 
ranges. The former term is not localized geographically, 
so that —kw, must represent the effect which is responsible 
for the lee-slope development. It is conceivable that the 
orographic contribution can influence the thermal pattern 
in such a way that —V,-V(ém-+-J) itself is altered, making 
that a secondary effect of the orography. Since , is 
positive for winds blowing from higher to lower levels, 
downslope winds will always contribute negatively to 
0¢/Odt and increase thermal ridging. 

The entire expression (22) is quite consistent with other 
types of manipulation of the two-parameter model at- 
mosphere, such as those employed by Sawyer and Bushby 
[8] and Sutcliffe [9], but is more appropriate for the present 
problem in that it refers explicitly to the thermal pattern. 
Equation (21) is not limited to orographic effects, and 
can be utilized to study any kind of large-scale thermal 
pattern. Consider a typical mid-latitude trough-ridge 
pattern at the mean level, where a low-level baroclinic 
cyclone appears on the eastern side of the upper trough. 
The with the surface Low will 
advance northward, becoming greatly extended, and 
finally become concentrically closed. Because of the 
close association between changes in the potential thermal 
vorticity and changes in the thickness, the first term on 
the right side of (21) explains the northward movement 
of the thermal pattern. Now, the origin of closed iso- 
therms cannot be explained by a single borizontal advec- 
tion process, and therefore, a second effect, such as that 
described by the last term on the right side of (23), must 
be responsible for the development of pools of cold and 
warm air. For this reason equation (21) could be appro- 
priate in studying cold Lows and warm Highs. 


isotherms associated 


3. THE IDEALIZED MODEL 


Consider an ideal sinusoidal zp field as schematically 
represented in figure 3. The z-axis lies in the east-west 
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Figure 3.—Schematic model of ideal thickness pattern. 
of Zmo in the z, y plane. z’’=axis of ridge. 2’= 
maximum descent. 


direction, and the y-axis coincides with the north-south 

mountain ridge. The crest of the idealized thermal ridge 

coincides with z’’. The following statement can then be 
made 

. Q2rz 

Zmo=a@ sin — +by, (24) 

where a is the amplitude of the ridge, \ the wavelength 

of the pattern, and 4 is taken to be a constant, although it 


generally slightly increases toward the north. Now if the 
geostrophic assumption holds, 


g g(2e\?_ . 2nz 
Smo=> V? z= —9 (=) asin =~: 


f Sf 


The above expression indicates the resemblance between 
@ and Zmo since 


(25) 


2 
o=a’ sin << +b'y, (26a) 


where 


(26b) 


a’ =— 8(7)+K |a, and b’=—K b. 


Consider a pattern of flow at the mean level, where the 
component of the mean wind u,,= U, and v,, is represented 
by a perturbation v’, and the term, —Vyno-V(fém+/) 38 
represented by the symbol Q. The 500-mb. pattern 
consists of a steady zonal current with a speed U’ upon 
which is superimposed a transversal wave perturbation. 
The expression (21) becomes 
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(27) 


In order to determine the eastward movement of the 
thermal ridge, this equation is differentiated with respect 
to x at the point 2x’’ 


2 (5).--Y [Fy UG) +K } Ja 


pm ov’ og 
he Ox x + 


(28) 
The last two terms on the right side of (28) represent the 
influence on the motion of the thermal ridge by the per- 
turbed part of the mean flow, and therefore demonstrate 
the effects of migratory troughs and ridges at the mean 
level. It is evident that Ov’/Or is merely ¢,,, and that the 
entire term, Ab¢,,, represents an advective influence by 
the non-zonal component of the mean wind. The last 
term represents a generating effect on thermal ridges by 
vorticity advection gradients at the mean level. Hence- 
forth, these two terms will be combined as one, where 


. , ot 
(’=Kb¢e,+ 


— From (28) the following conditions exist : 


7 (ar) 


~) =0, stationarity 2 
dr Var ),-. stationarity (29a) 


° (2*) (29b) 
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First consider the equation of stationarity: 
‘’, (28) becomes 


—*(& rte 


1 “oly ) {47 (=) y+x}] 


If the magnitude of Q’ is less than the magnitude of 


Ow, Ow, ‘ ol 
k ey “must be a negative quantity in order for a to be 


positive, 


>0, quasi stationarity 
- 


<0, eastward movement (29¢) 


solving for a 
at sr 


Thus stationary thermal ridges must become 


positioned where = is negative, in particular at a point 


downstream from the region of maximum descent. 
Equation (30) states that there exists a positive finite 
amplitude which will satisfy the stationarity condition. 
This amplitude will be called the “critical amplitude” 
represented by the symbol a,. Substitution of (30) in 


(28) vields 
g 2a P , — 
§(57) +x } la. a)- (31) 


ate" L(x) 


or (3 
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Equation (31) explicitly states that the thermal ridge will 
advance eastward, if a,>a. 

Consider the idealized development of a thermal ridge 
when a steady pattern of w, is imposed upon the lee slope 
of the mountain at =0. The initial position of the ridge 
will be at z’, since the generating effect of —kw, is greatest 
at that point. Now 0w,/dz is zero at 2’, and therefore, 
according to (30) and (31), the ridge must move east- 
ward in the first instant, but as it does so, a, also increases 
according to (30). The term, —kw,, represents an effect 
which is continuously generating a thermal ridge, and 
thus increasing the amplitude of the ridge with time. 
The critical amplitude can never exceed the actual am- 
plitude without imposing the stationarity condition. 
Similarly, the actual amplitude can never greatly exceed 
the critical amplitude, since that quantity will also in- 
crease as the thermal ridge moves east of zx’. If a=a,, 
equation (31) states that an approximate condition of 
stationarity exists. Although the extent and duration of 
approximate stationarity cannot be shown, the condition 
will persist until the thermal ridge attains a position 
downstream from 2’, where 0w,/dr is at a minimum. 
Further eastward movement will be accompanied by a 
rapid decrease in critical amplitude according to (30), and 
a negation of the condition of approximate stationarity. 

The horizontal advection of ¢ represents a decrease in 
this quantity ahead of z’’, and an increase behind it, while 
—kw, is decreasing ¢ ahead of z’’, but is decreasing it more 
strongly at z’. The growing importance of the advection 
of ¢ causes the ridge to move slowly eastward. The ulti- 
mate reduction of the critical amplitude corresponds physi- 
cally to the domination of the horizontal advection over 
the orographic process. 

If Q’ is negative across the ridge axis, the numerator in 
(30) will tend toward zero, leaving aa. According 
to (31), this event will tend to remove an approximate 
stationarity condition which may be present. If no ridge 
is present, none may be allowed to form under such cir- 
cumstances, since any amplitude which is created will be 
immediately advected downstream. Similarly, when Q’ 
is positive, equation (30) states that the maximum critical 
amplitude will be much larger than in the situation where 
(’ is negative. In effect, this means that the approximate 
stationarity condition exists for a considerably longer 
time. 

Since 6 is a negative constant in the present reference 
system, from the definition of Q’, negative relative vortic- 
ity at the mean level contributes to stationarity. The 
generative effect of Q will also contribute to stationarity 
when 0Q/Oz is positive, or when stronger cyclonic vorticity 
advection occurs east of the thermal ridge. Calculations 
of this term have shown its magnitude to be very small, 
and therefore Q’ primarily consists of [,. Thus Q’ will 
tend to be strongly positive in the vicinity of a 500-mb. 
ridge, and begin to decrease sharply behind it. 

In conclusion, it must be pointed out that (1) the initial 
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1000-mb. height contours (ft.), 1200 Gmv, January 12, 
1959 (unit digits omitted). 


Figure 4. 


formation of a thermal ridge is dependent on the funda- 
mental assumption that appreciable downslope velocities 
exist, and (2) the initial formation of the thermal ridge 
will not be observed as a lee-slope phenomenon if large 
positive values of relative vorticity at 500 mb. are occur- 
ring over the lee slopes of the mountain. If Q’ is the 
largest in the numerator of (30) the position of the thermal 
ridge will not necessarily be localized on the east slope of 
the mountain range. Finally (3), the magnitude of both 
terms in the numerator of (30) will tend to vary with the 
wind speed, and therefore the critical amplitude may be 
somewhat insensitive to wind speed. 


4. EXAMPLES OF LEE CYCLOGENESIS 
A SPECIFIC CASE 


An example of pseudo front development on the lee 
slopes of the Rocky Mountains has been chosen. The 
situation selected occurred between 1200 Gut, January 12, 
1959, and 1200 emt, January 13, 1959. The series of 
surface maps (figs. 4, 7, 10) reveals the development of 
the surface front. It was first noticed in connection with 
a feeble trough in the 065 contour in figure 4, and later 
exhibited a pronounced wind shift line, although no 
marked horizontal temperature contrast’ was present. 
The series of thickness maps (figs. 5, 8, 11) indicates the 
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Figure 5.—1000-500-mb. thickness contours (ft.), 1200 Gt, 
January 12, 1959 (ten of thousands and unit digits omitted). 
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500-mb. height contours (ft.), revealing strong eyclonic 
vorticity advection from the Pacific, 1200 Gm7, January 12, 1°59 
(ten of thousands and unit digits omitted). 
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Figure 7.—-1000-mb. height contours (ft.), 0000 Gmv, January 13, 
1959 (unit digits omitted). 


typical sinusoidal shape of the thermal pattern. Of note 
is the deep pocket of warm air around Colorado, associated 
with strong descending motion in that area. Although 
not presented here, maps preceding 1200 emt, January 12, 
indicate the presence of a weak thermal ridge. Over 
Canada, the temperature field associated with the vigorous 
Pacific Low has apparently merged with the northern 
portion of the thermal ridge, and continued to move 
rapidly eastward. In terms of equation (30), the large 
evelonie vorticity advection at 500 mb. (fig. 6) over south- 
ern Canada was removing the stationarity condition over 
the eastern slope of the mountain. South of the Canadian 
border, the 500-mb. ridge was positioned slightly west of 
the thermal ridge (figs. 9, 12). As the 500-mb. ridge 
moved east of the thermal ridge, the latter began to move 
eastward. This event could be anticipated by the idealized 
model, since Q’ would tend to decrease sharply upon the 
passage of a 500-mb. ridge. Large Pacifie Ocean cyclones 
Which approach Puget Sound create strong westerly flow 
across a considerable part of the mountain chain. Ac- 
cording to the model, this event occurring simultaneously 
with the approach of a 500-mb. ridge from the west, is 
ideal for the establishment of a stationary thermal ridge 
on the lee slopes. 
GENERAL CASE 


A general approach was utilized to study lee cyclogene- 


sis. All instances of thermal ridge intensification on the 
lee slopes of the Rocky Mountains were systematically 














Figure 8.—1000-500-mb. thickness contours (ft.), OOOO Gm, 
January 13, 1959 (ten of thousands and unit digits omitted). 
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500-mb. height contours (ft.), 0000 cmv, January 13, 
1959 (ten of thousands and unit digits omitted). 
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Figure 10.—1000-mb. height contours (ft.), 1200 Gm7, January 
13, 1959. A-~B is position of atmospheric cross-section repre- 
sented in figure 2 (unit digits omitted). 


tabulated over the months December through March 
1958-59 and 1959-60. The classification of the ridges 
stationary” and ‘‘mov- 
ing’. The decision concerning the category of motion 
was qualitative, and based upon whether the thermal 
ridge appreciably moved eastward during the following 
12-hour map interval.?. Thermal ridges which remained 
stationary on more than one successive map were tabu- 
lated as “stationary”? more than once. If the ridge later 
began to move eastward, it was tabulated again as ‘“‘mov- 
ing.”’ For each case the position of the 500-mb. ridge 
and trough, relative to the thermal ridge, was tabulated. 
The results of this undertaking are presented in figure 13. 
The outstanding feature is the shift of the mean position 
of the 500-mb. ridge in the ‘‘stationary’”’ cases to the east 
of the thermal ridge in the ‘‘moving” ones. Similarly, 
the mean position of the 500-mb. trough is situated at a 
distance greater than 15° of latitude west of the “station- 
ary”’ thermal ridges, and about 12° west of the “moving” 
ones. The mean direction of the 500-mb. wind across the 
axis of the thermal ridge was 290° for the ‘‘stationary”’ 
” thermal ridges. In 
all of the cases of thermal ridge development and motion 


oe 


was divided into two categories, 


ae 


cases and 265° in cases of “moving 


(about 55), only four instances were accompanied by no 
definite trough-ridge pattern at 500 mb. The importance 


2 The surface map was always categorized before the 500-mb. flow was observed and 
categorized. 














FIGURE 


11.—1000-500-mb. thickness contours (ft.), 1200 Gaz, 


January 13, 1959 (ten of thousands and unit digits omitted). 
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FIGURE 


12.—500-mb. height contours (ft.), 1200 emr, January 13, 
1959 (ten of thousands and unit digits omitted). 
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Figure 13.—Occurrence of 500-mb. troughs and ridges in 5° lati- 

tude intervals with respect to the thermal ridge axis, for station- 
Compiled for the winters of 1958-59 and 
1959-60 (units of ordinate go from 0 to 12 cases). 


ary and moving cases. 


of the 500-mb. ridge in influencing the motion of lee-side 
disturbances is expressed by the term Q’ in the numerator 
of equation (31). As the 500-mb. ridge approaches the 
lee slopes of the mountains, not only are conditions favor- 
able for the initial development of thermal ridges, but 
those ridges which form will remain stationary for a short 
time due to the increasing effect of Q’ which offsets the 
tendency for the thermal ridge to be advected downstream. 
Figure 13 also suggests that in the vicinity of 500-mb. 
troughs the numerator of (31) becomes vanishingly small, 
preventing the formation of thermal ridges. 


5. SUMMARY AND CONCLUSIONS 


The ultimate purpose of this paper was to study the 
development and motion of the frequently observed lower 
tropospheric frontogenesis associated with thermal ridge 
intensification on the lee slopes of the Rocky Mountains. 
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The dynamical formulation showed that changes in the 
thermal pattern can be described by a potential thermal 
vorticity equation, which consists of three terms: (1) a 
term representing the advection of the potential thermal 
vorticity by the mean wind; (2) an orographic term; and 
(3) a term representing the advection of absolute vorticity 
at the mean level by the thermal wind. The mean level 
was considered to be 500 mb. An idealized sinusoidal 
model of the thermal ridge was devised in order to demon- 
strate the effect of each of the three terms on its motion. 
The results showed that the orographic effect was likely 
responsible for the thermal ridge amplification on the lee 
slopes of the mountains. According to the idealized 
pattern, the thermal ridge remains approximately sta- 
tionary downstream from the point of maximum descent 
during the intensification process. Its amplitude con- 
tinues to intensify until it attains a maximum critical value. 
At this point, the advection of potential thermal vorticity 
by the zonal flow at 500 mb. begins to dominate the mo- 
tion, and the thermal ridge accelerates eastward. The 
effects of the perturbation in the 500-mb. current were 
shown to highly influence the development and motion of 
thermal ridges. The presence of migratory troughs and 
ridges at 400 mb. either (1) increases the maximum critical 
amplitude for stationarity, and thus increases the length 
of time for development, or (2) decreases the maximum 
critical amplitude, and therefore removes the stationarity 
condition, or (3) prevents any localized lee-slope thermal 
amplification from occurring. 

An actual case was chosen to examine the development 
and motion of a thermal ridge in terms of the dynamical 
considerations. Th's undertaking provided a link be- 
tween the theory and actual synoptic experience. 

A survey of many cases was taken over the period of two 
winters. The results presented support the theory. 
Almost every case selected, whether it was stationary or 
moving eastward, was associated with a ridge at 500 mb. 
In the moving cases the 500-mb. ridge was generally 
centered east of the thermal ridge, and in the stationary 
instances it was centered west of the thermal ridge. 
Combined with the theory, these results lead to some 
conclusions concerning the nature of pseudo front develop- 
ment and motion. 

The forecaster should consider the following criteria: 
(1) The thermal ridge will develop when the surface flow 
is such as to produce large-scale descent on the lee slopes 
of the mountains. (2) No thermal ridging will take place 
when the 500-mb. ridge lies downstream from lee slopes. 
(3) Thermal ridging will appear with the approach of a 
500-mb. ridge from the west, and (4) will move eastward 
upon the passage of the 500-mb. ridge. 
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ABSTRACT 


Solar and terrestrial radiation measurements that were obtained at Amundsen-Scott (South Pole) Station and 


on Ice Island (Bravo) T-3 are presented for representative summer and winter months. 


Of the South Polar net 


radiation loss during April 1958, approximately 20 percent of the energy came from the snow and 80 percent from the 


air. 
rate. 


The actual atmospheric cooling rate during that period was only about 1/6 of the suggested radiative cooling 
The annual net radiation at various places in “Antarctica is presented. 


During 1958, the South Polar atmos- 


phere transmitted about 73 percent of the annual extraterrestrial radiation, while at T-3 the Arctic atmosphere 


transmitted about 56 percent. 


The albedo of melting sea ice is discussed. 
indicate that the net radiation is positive on both clear and overcast days but greatest on overcast days. 


Measurements on T-3 during July 1958 
Refreezing 


of the surface with clear skies, as observed by Untersteiner and Badgley, is discussed. 


1. INTRODUCTION 


The elliptical orbit of the earth brings it about 3 million 
miles farther from the sun at aphelion than at perihelion; 
consequently, during midsummer, about 7 percent less 
solar radiation impinges on the top of the Arctic atmos- 
phere than on the Antarctic atmosphere during a com- 
parable period. This difference is enhanced as _ solar 
energy penetrates into both polar atmospheres. Absorp- 
tion, scattering, and reflection of the solar rays gives each 
polar region its own particular radiation environment. 
There are also other notable differences between the heat 
budgets of these two areas; for example, the conduction of 
heat through the ice is distinctly different. Annually, 
heat from the Arctic Ocean is conducted upward through 
the thin ice pack to the relatively cold surface where the 
temperature averages about —20° C. In contrast, the 
flux of heat through the ice layers of central Antarctica 
isquite small. Because of the heat budget differences, the 
annual temperature near the North Pole is about 30° C. 
warmer than that at the South Pole. It is the purpose of 
this paper to discuss some aspects of the thermal energy 
budgets of these regions. 

The data which are presented were obtained during the 
International Geophysical Year and later years at Amund- 
sen-Scott Station, located within a mile of the geographic 
South Pole, and at Ice Island T-—3, drifting in the Arctic 


_ 


‘Paper presented at the International Antarctic Symposium at Buenos Aires, No- 
vember 1959, 


Ocean. This Island was about 5 by 11 miles in size and 
about 52 meters thick (Crary et al. [2]) in 1953 when it 
drifted near 88° N., 100° W. In the years that followed, 
this it drifted southward and in July 1958 was located 
79.5° N., 118° W. 

Solar radiation measurements at both stations were 
obtamed. with Eppley pyranometers. The data are 
corrected for the temperature response of the instrument 
(MacDonald [11]) and are presented in the International 
Pyrheliometric Scale of 1956. A‘ both stations, Beckman 
and Whitley (Gier and Dunkle type) radiometers were 
used to measure the combined solar and_ terrestrial 
radiation streams. 


2. WINTER MONTH AT THE SOUTH POLE 


With the exception of a few weeks of twilight, sunset 
at the time of the March equinox marks the beginning of 
6 months of continuous darkness at the South Pole. 
During the first of the dark months, April, the tempera- 
ture a few meters above the snow averaged —58° C. 
(1957—59)—the same as the average temperature during 
the entire dark period. 

In April 1958 the long-wave radiation ? from the snow 
surface averaged 229 ly. day~' (table 1), while the atmos- 
pheric (back) radiation returned 76 percent of this energy 
(175 ly. day") to the surface. The net radiation aver- 
aged —54 ly. day. 


? Terrestrial radiation, with most of the energy between wavelengths of 3 to 30 microns. 











Thermal energy exchange at the snow surface, Amundsen- 
Scott (South Pole) Station, Antarctica 


Annual 
1958 


January April 
1958 1958 


Incident solar radiation (ly. day~')* : . 770 0 266 
Reflected solar radiation (dy. day~')_. 674 0 218 
Albedo (percent) xs ™ 82 
Snow surface radiation (ly. day~) ; 440 229 301 
Atmospheric (back) radiation (ly. day~')-. 309 175 218 
Return (percent) 70 76 72 
Net radiation (ly. day) —35 —i4 —35 
Thermal energy from snow (ly. day~!) cece 11 | act 
Thermal energy from air (ly. day~') aS Fee 43 |. pio 


*One Langley (ly.) equals one cal. em.~? 






The total thermal energy in the top 12 meters of snow 
at the beginning and end of April 1958 was calculated, 


using the equation 


"12 
Q= } cp Tdz (1) 


0 


where ¢ is the specific heat in cal. gm~! deg.~' (List [9}), 


p is the density in gm. cm.~* (Giovinetto [4]), and 7 


is the temperature (fig. 1) in °K. The calculations 
indicate that during the month heat was conducted to 
the surface at the rate of 11 ly. day~'. This suggests 
that, of the net radiation loss during April, about 20 per- 
cent of the energy came from the snow and 80 percent 
from the air. During clear, cold periods at the South 
Pole the snow and air supply about equal amounts of en- 
ergy to make up the surface radiation loss (Hanson [6]). 
Liljequist [8] found that with clear skies at Maudheim, 
along the coast of Antarctica, roughly 40 percent of the 
required energy comes from the snow and 60 percent from 
the air. 

Some idea of the heat budget of the atmosphere during 
this period can be obtained from the airborne radiation 
measurements which taken with Suomi airborne 
radiometers (Suomi et al. {[14]). Data from the clear-sky 
flight on April 27, 1959 (fig. 2) indicate a radiative loss of 
240 ly. day! at 50 mb. Assuming this loss is representa- 
tive of April 1958, and adding the small amount of heat 
which was conducted from the snow (11 ly. day~'), the 
net cooling rate from the surface to 50 mb. becomes 
1.49° C. day~'. This is about 6 times greater than the ob- 
served cooling rate (0.26°C. day~'). Presumably, subsi- 
dence and advection provide the necessary energy to ac- 
count for the discrepancy. 


were 


3. SUMMER MONTH AT THE SOUTH POLE 


Even though the South Polar plateau receives more 
solar radiation at midsummer than any other area on 
earth, the temperature of the snow surface remains well 
below freezing. In January, the warmest month of the 
summer, the temperature averages near —27° C., and 
rarely exceeds about —17° C. 

During January 1958, with continuous sunlight at the 
South Pole, the incoming solar (sun and sky) radiation 
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Figure 1|.—Snow temperature profiles of April 1 and 30, 1958, and 
monthly temperature-change profile, Amundsen-Scott (South 
Pole) Station, Antarctica. 




















averaged 770 ly. day~'. Measurements indicate that, of 
this amount, SS percent (674 ly. day~') was reflected from 
the snow, while the remainder (96 ly. day~!) was absorbed. 
During the same period, the snow surface radiated 440 ly, 
day~', while the atmospheric (back) radiation returned 70 
percent (309 ly. day~') of that energy. The net radiation 
averaged —35 ly. day™'. 


4. ANNUAL ENERGY EXCHANGE, ANTARCTICA 


During the 6 months of sunlight at the South Pole, the 
solar radiation which was incident on the snow totaled 
9.7110* ly. The snow reflected 7.96 104 ly., indicat- 
ing an average albedo of 82 percent for the sunlit period. 

Unlike solar radiation, the emission of long-wave radia- 
tion by the snow is continuous throughout the year. Dur- 
ing 1958, the snow surface radiation averaged 301 ly. day~, 
of which about 72 percent (218 ly. day~') was returned by 
atmospheric (back) radiation. This percentage is rela- 
tively unchanged from summer to winter even though the 
sky is much clearer during winter (fig. 3). With other 
things being equal, clear skies would certainly tend to lower 
this percentage. Apparently, a compensating factor is 
that the surface temperature inversion is more intense 
during winter; this would allow a greater return of the 
surface radiation. 

The net radiation at the South Pole (2800 m.) averaged 
about —35 ly. day~! during 1958. Liljequist [7] found an 
annual loss of 25 ly. day“! at Maudheim, and Loewe [9 
found a loss of 20 ly. day~! at Port Martin. Both stations 
are located along the coast of Antarctica. Rusin [13] has 
reported an annual net radiation of —6 to —8 ly. day‘ at 
Mirny, another coastal station, and —19 to —22 ly. day™ 
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Figure 2.—Net radiation profile, temperature profile, and tempera- 
ture-change profile computed from radiometersonde (Suomi 
et al. [14]) ascent at Amundsen-Scott (South Pole) Station, 
Antarctica, 0600 Gur April 27, 1959. 


at Pionerskaya (2700 m.) in the interior of eastern Ant- 
arctica. There is little doubt that, as a whole, the surface 
of Antarctica has a negative net radiation, although cer- 
tain snow-free areas on the continent and portions of the 
Palmer Peninsula are probably positive. Considering the 
magnitude of the individual incoming and outgoing radia- 
tion streams, the slight variation of net radiation, as ob- 
served in Antarctica, seems rather remarkable. 


. ANNUAL ENERGY EXCHANGE, ARCTIC OCEAN 
BASIN 


Although the elliptical orbit of the earth causes com- 
paratively less midsummer extraterrestrial radiation in 
the Aretic, it also provides 9 additional days of sunlight 
at the North Pole each year compared to the South Pole. 
The net result is that during the course of their respective 
sunlit periods, equal amounts of solar energy impinge on 
the top of the North and South Polar atmospheres. This 
Was pointed out by Milankovitch [12], who indicated that 
at both geographic poles the annual extraterrestrial radia- 
tion is of equal intensity, and, assuming a solar constant 
of 2.00 ly. min.~!, totals 13.3310 ly. during the year. 
This equality is not maintained as the solar rays penetrate 
the polar atmospheres, however. The previously men- 
tioned measurements at the South Pole indicate that, 
annually, about 73 percent of the extraterrestrial radiation 
Was incident on the snow surface. At Ice Island T-3, on 
the other hand, only about 56 percent of the extraterres- 
trial radiation * was incident at the surface during 1958. 


_—— 


'Lo (dat 79.5° N., the extraterrestrial radiation at T-3 totaled 13.81X 10 ¢ ly. during 
1958 
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Figure 3.—Average cloud coverage during the period January 1957 
through December 1960 at Amundsen-Scott (South Pole) 
Station, Antarctica. 


The fact that a smaller percentage of the extraterrestrial 
radiation reaches the surface in the Arctic may be due to 
a number of factors: possibly there are thicker cloud sys- 
tems in the Arctic, the surface albedo is less in the Arctic, 
and the optical thickness of the Arctic atmosphere is 
greater because of the comparatively lower surface eleva- 
tion of the Arctic. Precisely how these variables affect 
the amount of solar energy, incident on the Arctic and 
Antarctic, is difficult to determine without a better knowl- 
edge of the surface albedo and thickness of the cloud 
systems (Fritz [3]). 


6. SUMMER MONTHS IN THE ARCTIC 


The snow that covers much of the Arctic sea-ice in 
early spring is gradually melted during May and June. 
As a result, pools of melt water form on the ice floes in 
late June and remain during July and sometimes August. 
These pools aid in lowering the surface albedo during this 
midsummer ablation period. Sverdrup [15] has indicated 
that the albedo of melting Arctic sea-ice is between 60 and 
65 percent. Recent Soviet investigations (Briazgin [1]) 
have shown a similar albedo value, 60 percent. The fact 
that pools of melt water on the floes lower the albedo is 
indicated by the results of an aerial albedo survey over 
the Arctic Ocean. This survey indicated that the albedo 
of melting ice with a maximum amount of puddling is 
about 46 percent (Hanson [5]). 

Measurements on T-3 indicate that during July 1958 
the incident solar radiation averaged 524 ly. day™'. 
Assuming an albedo of 60 percent, the ice would have 
gained 210 ly. day~'. Measurements also indicate that, 
with an average 6.7-tenths cloud coverage, the atmospheric 
long-wave radiation returned about 88 percent (578 ly. 
day~') of the 653 ly. day~' which were emitted as long- 
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wave radiation from the surface. This gives an average 
net radiation of +135 ly. day~ during July. 

Currently, one of the most important questions in polar 
heat-budget investigations is: How is the “surplus” 
thermal energy, available from radiative exchange, used 
in warming or melting the ice, evaporation, or warming 
the lower atmosphere? During the previously mentioned 
period, for example, the net radiation could have melted 
58 em. of surface ice, assuming a density of 0.9 gm. em.~* 
and a latent heat of fusion of 80 cal. gm.~' The actual 
ablation was probably somewhat less, however, as a small 
amount of energy is lost in evaporation (Untersteiner and 
Badgley [16]), and possibly, as Fritz [3] and Yakovlev [17] 
have suggested, some energy may be lost to the atmosphere 
and 
temperature profile measurements are not available, the 
heat and water budget during this important summer 
ablation period cannot be determined precisely. 


by turbulence. Because ablation, evaporation, 


An interesting observation regarding the heat budget 
was made by Untersteiner and Badgley [16] on Floating 
Ice Station “A”: 


During the summer, melting oecurred mostly when there was 
overcast and strong atmospheric radiation. Radiosonde observa- 
tions reveal frequent inversions, with comparatively high cloud 
temperatures. Temporary breaks in the overcast were frequently 
accompanied by freezing at the surface even though direct solar 


radiation was relatively larger in such periods. 

In order to investigate the role of radiative heat in melting, 
we have examined the radiation data for clear and over- 
cast conditions at T-3 during July 1958. 
ments of the incoming solar and atmospheric (back) radi- 


From measure- 


ation, assuming an albedo of 60 percent which is probably 
representative of the melting pack ice, it was found (table 
2) that on 6 clear days the net radiation averaged +76 
ly. day~'. On 9 overcast days it averaged +147 ly. 
‘an increase by a factor of about 2 providing the 
albedo is unchanged. The radiation data, as presented 
here, lend some support in explaining these observations 
made by Untersteiner and Badgley [16]. However, a full 


day 


TABLE 2.— Thermal energy exchange at the surface of Arctic sea-ice 
during clear, overcast, and average sky conditions as determined 
from measurements of incident solar and atmospheric (back) radia- 
tion at T-3, assuming an albedo of 60 percent and ice surface 
temperature of O° C. 


July 1958 


Clear sky 
condition, 
| <I-tenth 
sky cover 
(on 6 days) 


Averace 
sky condi- 
tion, 6.7 

tenths 
sky cover 
(31 days) 


Overcast 
sky condi- 
tion, 10- 
tenths 
sky cover 
(on 9 days) 


Incident solar radiation (ly, day ie 401 | 
Reflected solar radiation (assuming an albedo of 

60 percent) (ly. day~!) , 386 241 
Ice surface radiation (ly. day! oa as 653 
Atmospheric (back) radiation (ly. day~!)_. ry O40 
Return (percent as 
Net radiation (ly. day +147 
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explanation would require some idea of the rates of evap- 
oration and turbulent heat exchange in the boundary 
layer. 

There is also an additional point of interest. In view 
of the observed refreezing of the surface with clear skies, 
it is rather surprising that the net radiation should be 
positive. An examination of the data revealed that an 
albedo of 72 percent or higher would be required in order 
to have a negative net radiation with clear skies. This 
seems excessive. It seems reasonable that if freezing of 
the surface occurs with a positive radiation balance, either 
one of two things is happening. Either the total heat 
budget is negative because of losses by evaporation and 
turbulence and therefore refreezing of the surface occurs; 
or, the total heat budget is positive and the “surplus” 
energy is used in melting just below the surface while 


freezing continues on the surface. The latter may be 


possible as the ice is partially transparent to solar radia- 


tion and opaque to the longer wavelengths of terrestrial 
origin. In order to present a realistic model of the proc- 


esses involved, additional field studies are desirable. 
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THE WEATHER AND CIRCULATION OF FEBRUARY 1961’ 


An Example of Attentuation in the Long-Wave Pattern 


L. P. STARK 


Extended Fore 


1. CONTRASTING REGIMES IN WINTER 1960-61 


Attenuation of the long-wave features in the western 
portion of the Northern Hemisphere in February 1961 
was in marked contrast to the patterns that had prevailed 
in the preceding winter months. In December [1] and 
January [2] pronounced amplification of the monthly 
mean 700-mb. flow pattern dominated the circulation and 
effectively controlled the weather in North America. 
Widespread cold conditions in December were repeated 
in January in the East and South in response to the strong 
ridge in western North America that extended from the 
Great Basin to the Arctic Ocean. Substantial changes 
in the circulation in February 1961 (fig. 1) produced 
strong damping of the long waves in the mid-troposphere. 
As a result the country was subjected to mild Pacifie air 
masses instead of cold continental air; the westerlies at 
mid-latitudes over North America were strengthened; 
and temperatures in the East became warmer than normal 
near the end of the winter season. 


2. SOME ASPECTS OF THE MONTHLY MEAN 
CIRCULATION 


The changes in circulation from January to February 
are best examined by referring to figure 2, the difference 
in 700-mb. height departure from normal between the 
two months. The greatest changes were in the Pacific, 
where most of the area underwent height rises in Febru- 
ary, and the deep, cyclonic center of January at middle 
latitudes weakened and moved northward. Consequently, 
the ridge in western North America could no longer be 
sustained in such strength as was noted in December and 
January. Heights fell throughout the area from January 
to February, with maximum change of —420 feet in the 
Yukon (fig. 2). 

The trough in eastern North America advanced to the 
central Atlantic as flattening of the pattern progressed 

Height anomalies increased over most 
North America (maximum of +310 ft. in 
curvature changed from 


downstream. 
of eastern 
New York) where contour 


1 The weather of March, April, and May 1961 will be described in the June, July, and 
August 1961 issues of the Review, respectively. 


ist Section, U.S. Weather Bureau, Washington, D.C. 


strongly cyclonic in January to slightly cyclonic in 
February. The result of these changes in and near North 
America was a decided increase in the average speed of 
the westerlies at mid-latitudes. The hemispheric zonal 
index (5° W.-175° E.) was 2 m.p.s. greater in February 
than it was in January. This increase occurred despite 
a regional decrease in westerlies from the central Atlantic 
to Europe where the meridional component of the flow 
dominated. 

In western Europe changes in the average circulation 
at 700 mb. were large. Diffluent, cyclonic flow in Janu- 
ary was replaced in February by a strong ridge (fig. 1) 
in which there were maximum height rises of 390 feet in 
Spain (fig. 2). The evolution of this ridge was simul- 
taneous with damping of the pattern in the eastern 
Pacific and western North America; its growth deserves 
further comment, based on 5-day mean 700-mb. charts 
one week apart (not shown). 

On the 5-day mean chart for February 2-6, a_ high 
center was located just to the west of Morocco with a 
flat ridge extending northward to Scotland. One week 
later the High was centered over Gibraltar and the ridge 
was still rather flat to the north, but the trough southward 
from Greenland had become quite deep near 40° NX. 
Meanwhile a positive height anomaly center of 390 feet 
had moved from Gibraltar to northern Spain and in- 
creased to 590 feet. 

During the week of February 16-20 the anomaly center 
reached a maximum of 760 feet above normal in France 
and was associated with a dominant blocking High. In 
the last week of February the High was observed over the 
Baltic Sea, still strong, with an omega-configuration well 
defined by a strong trough (negatively tilted) in the east- 
ern Atlantic and a deep resonant trough from northwestern 
Siberia to the Red Sea. 

Another blocking High was observed on the 30-day 
mean chart (fig. 1) over northeastern Siberia, but this 
was considerably weaker than the one in Europe. Its 
development occurred as the trough along the Asian coast 
deepened at middle and low latitudes and as the ridge 
which was in the Arctic Ocean in January retrograded 
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Figure 1, 


Mean 700-mb. contours (solid) and departures from normal (dotted), both in tens of feet, for February 1961. 


The trough 


in the Southwest intensified as the wavelength and westerlies increased at mid-latitudes. 


Ordinarily the effect of blocking shows up in a migration 
southward of the westerlies and a decrease in the average 
zona! wind speed. But the blocking Highs this month 
were too far removed to have an appreciable effect on 
the zonal index in the western portion of the Northern 
Hemisphere or on the weather and circulation in North 
America. Dominant in the weather of the United States 
were the generally westerly flow from the eastern Pacific 
to the western Atlantic, and especially, the position of the 
Aleutian Low at an extreme northeast location. 


3. WEATHER RELATED TO THE MEAN CIRCULATION 


TEMPERATURE 


Departures of average surface temperature from normal 
(fig. 3) were predominantly positive except for a small 


area in the Southwest where it was 1°-2° F. cooler than 
normal. The correspondence was close between the flow 
pattern and the temperature anomaly. Compare the 
height departure from normal (dotted lines fig. 1) with 


the temperature departures in figure 3. Considering 
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Figure 2.—Change of mean 700-mb. height departures from normal 
(tens of feet) from January 1961 to February 1961. Attenua- 
tion of flow pattern near North America in February is accented 
by falls in the ridge in western North America and rises in the 
troughs in the eastern Pacific and western Atlantic. 
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Departure of average temperature from normal (°F.) 
Height anomaly (fig. 1) and temperature 
(From [5].) 


Figure 3. 
for February 1961. 
anomaly are almost identical in pattern. 


sign only, the similarity in the patterns is remarkable; 
considering magnitude, the analogy is weaker. 

The correlation between local 700-mb. heights and 
surface temperature anomaly is generally positive, but 
average correlations in the winter months explain only 
about 35 percent of the temperature variance [3]. By 
statistically integrating the behavior of the atmosphere 
at various distances from the local area, the explained 
variance can be considerably increased [4]. The synoptic 
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TABLE 1.—Class changes of monthly average surface temperature 
anomalies in 100 cities from January to February 1961 





Frequency 
(percent) 


Class change* | 


*+indicates a change toward higher average temperature. Classes are much below 
normal, below normal, near normal, above normal, and much above normal, 


interpretation of this statistical result is dependent on 
the well-known fact that the position of the Aleutian Low 
has a controlling effect on temperatures in much of North 
America. In February it is normally located over 
Kamchatka. When it is displaced far to the east (as it 
was in February 1961), 700-mb. heights are usually lower 
than normal in northwestern North America (as they 
were this month)(fig. 1). The resulting westerly flow 
prevents outbreaks of Arctic air from penetrating into 
latitudes below 40° or 50° N. Furthermore, increased 
westerlies at these latitudes in the eastern Pacific (4-8 
m.p.s. stronger than normal this February) undergo 
strong katabatic warming after striking the middle 
Rockies. The total effect is frequently unseasonable 
warmth east of the Divide at middle latitudes. 

The change in temperature regime from January te 
February was determined by the circulation change 
(fig. 2). Consequently, persistence of temperature was 
large west of the Mississippi River, where height changes 
were small; to the east, temperatures changed radically 
where height changes were relatively large. A summary 
of the January to February changes in monthly average 
temperature (in terms of five classes) at 100 cities (table 1) 
indicates the following: (a) 74 percent of the temperatures 
became warmer, (b) 7 percent became cooler, and (¢) 
21 percent became warmer by three classes (from much 
below normal to above normal). 


PRECIPITATION 


Precipitation in February (fig. 4) was much more 
extensive and generally heavier (in 64 of 100 cities) than 
it was in January [5]. From the Ohio Valley to the 
western Plateau and from the Canadian border to the 
Southern Plains, little or no precipitation fell in January, 
when northerly flow at 700 mb. and sea level inhibited 
the transport of moist air masses into the country. In 
February the average flow aloft (fig. 1) and at sea level [6) 
encouraged the advection of moisture northward from 
the Gulf of Mexico and eastward from the Pacific. 

In the Pacific Northwest there was a deep onshore flow 
from sea level to 700 mb. This current, the core o 
maximum wind speed or “jet stream’ at 700 mb., sup 
ported a very high frequency of frontal systems, 8 to 2 
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crossing the coast in February. As much as 20 in. of 
rain fell along the Washington coast as moist, Pacific air 
was forced to ascend the mountain slopes. Precipitation 
records for February were broken at Olympia, Wash. 
(7.02 in. above normal), Sexton Summit, Oreg. (6.75 in. 
above normal), and several other stations in Washington 
and Oregon. Flooding was reported in the Puget Sound 
area, the Olympic Peninsula, and the Willamette Basin. 

In the Southeast the distribution of precipitation was 
related almost ideally to the 700-mb. pattern [7]. Heavy 
precipitation occurred in the southwesterly flow ahead 
of the trough and through the weak ridge (fig. 1). There 
was additional support from the southeasterly anomalous 
component of the 700-mb. flow and from the warm, moist 
air at sea level that had a long trajectory from the Tropics. 
Locally severe flooding from Mississippi to Georgia 
resulted from 8 in. to more than 17 in. of precipitation that 
fell in a narrow band from Louisiana to western South 
Carolina (fig. 4). Several cities in Alabama and Georgia 
reported the wettest February of record. From early 
reports precipitation was heaviest in Birmingham and 
vicinity, where the total was 17.67 in., 12.50 in. greater 
than normal. 

In contrast to the wet weather just described were 
continued drought conditions over most of the Southwest. 
Many cities had near-record dryness for February and 
for the winter season. At Cheyenne, Wyo., and Grand 
Junction, Colo., February was the 16th and 11tb consecu- 
tive month, respectively, in which precipitation was less 
than normal. <A deficit of precipitation in February 
contributed to the driest winter of record at Chicago, IIl., 
Lansing, Mich., and Burbank, Calif. 


4. INTRA-MONTHLY TRANSITION 
JANUARY 31-FEBRUARY 14 


The 700-mb. circulation in the first half of February 
(fig. 54) was transitional between the large-amplitude 
circulation of January and the higher index state during 
the last half of February. 

In the eastern Pacific, despite appreciable filling, the 
eastern lobe of the Aleutian Low was still rather strong. 
The sea level counterpart (fig. 5B) was 11 mb. below 
normal, as intense as any vortex on the chart, but con- 
siderably weaker than the same Low in January, when 
pressures for the month averaged 20 mb. below normal [2]. 

The ridge over North America had much less amplitude 
than in January and was somewhat farther east. Height 
falls at 700 mb. eroded the western portion of the ridge at 
middle and high latitudes as the center of positive height 
departure from normal jumped from British Columbia in 
January to Hudson Bay in the first half of February. 

In the Atlantic the strong trough at 700 mb. progressed 
from the east coast and deepened as cold air was injected 
from Canada. The Icelandic Low at sea level was near its 
horm:| position, but was 11 mb. deeper than normal. 

Th» distribution of temperature departures from normal 
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Figure 4.—Total precipitation, inches, for February 1961. Locally 
severe flooding in the Northwest and in the Southeast was closely 
related to strong southwesterly flow in advance of mean troughs. 
(From [5].) 


for January 31—February 14 is shown in figure 5C. From 
the Central Plains westward, temperatures ranged from 
4°-18° F. warmer than normal. Most of this warmth 
was a function of the ridge aloft; some was related to the 
increasing westerlies and the deep Low in the Gulf of 
Alaska. In the eastern portions of the country cold air 
was closely associated with the anomalous northerly 
components of the flow at 700 mb. and at sea level. The 
source and trajectory of the cold air observed in the East 
was implicit in the strong ridge at sea level from Hudson 
Bay to the Gulf Coast (fig. 5B). 


FEBRUARY 15-MARCH 1 


During the last half of February broadscale damping 
of the major perturbations occurred in mid-latitudes 
from the central Pacific to the eastern Atlantic (fig. 6A). 
The Aleutian Low continued to weaken as heights in- 
creased as much as 480 feet in the eastern Pacific. At 
the same time heights decreased over Canada, and cy- 
clonic flow prevailed from Alberta to Labrador, in place 


of the ridge formerly there. In the western Atlantic 
height rises were widespread and reached a maximum of 
390 feet where there had been a trough the first half of 
the month. This nearly simultaneous growth of the 
Bermuda High and the High in the eastern Pacific suggests 
a teleconnection often noted in a period of rising index. 

General height falls in higher latitudes and rises in 
middle latitudes contributed to a pronounced strengthen- 
ing of the zonal index in temperate latitudes. At mid- 
month the 5-day mean 700-mb. zonal index (5° W.-175° 
E.) was 10.5 m.p.s., only 0.3 m.p.s. above normal. By 
the end of February, however, the index had jumped to 
13.6 m.p.s., 3.8 m.p.s. above normal. 

Deepening of the trough in the lee of the Rockies was 
assisted by the increased westerlies impinging on this 
natural barrier. Accompanying the increased cyclonic 
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flow were vigorous frontal systems that crossed the 
Rockies and intensified just to its lee. In fact, there 
were 16 days of frontal activity in Montana during the 
month, a maximum for the country. 

The sea level chart for February 15—March 1 (fig. 6B) 
should be compared with that of the first half of February 
(fig. 5B). Changes indicated are similar to those found 
on the 700-mb. charts for the same periods. The com- 
plementary anticyclogeneses in the eastern Pacific and 
the western Atlantic are clearly portrayed by the dotted 
lines showing the sea level pressure departures from 
normal. In the same period sea level pressures fell 5-9 
mb. from the Gulf of Mexico to Baffin Bay. The direc- 
tion of the flow and anomalous components are 
indicative of the changed regime. Note especially the 
northerly direction in the West and the southerly direction 
in the East. 


its 





Figure 5.—(A) Mean 700-mb. contours (solid) and departures from 
normal (dotted), both in of feet. (B) Mean sea level 
pressure (solid) and departures from normal (dotted), both in 
millibars. (C) Departure of average temperature from normal 
(°F.). All charts for January 31—February 14, 1961. 


tens 


A great warming in the East and a cooling in the West 
were the principal reactions to the reversed flow. Figure 
6C shows that average temperatures for February 15 
March 1 were generally 4°-8° F. warmer than normal 
over most of the country and a few degrees cooler in the 


Southwest. But the changes from the first half o! 
February to the last half were quite large. In the East, 
for example, temperature anomalies showed _ positive 


changes of 12°-16° F. in a belt from Florida to Nev 
York; west of the Plains States negative changes rangeé 
from 4°-10° F. 

As the transition in the circulation proceeded, ther 
was a definite change in the tracks of cyclones in and nea! 
North America. This might be anticipated from a stud} 
of figures 5 and 6, but for additional emphasis and clarity 
figure 7 is presented. This figure was prepared )y 
separating the tracks of migratory cyclones published ! 
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Figure 7.—Paths of migratory cyclones (A) Feb. 1-14, 1961, and (B) Feb. 15-28, 1961. Note general shift of cyclones northward 
both oceans and increased frequency of tracks originating in the Southwest from the first half to the last half of February. 
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[6] into halves; those originating in the first half of 
February (A) and those originating in the second half (B). 

These tracks clearly illustrate the changes discussed 
above. The transfer of cyclonic activity in the Pacific 
to a more northerly locus was a typical result of anti- 
cyclogenesis in the eastern Pacific. In the Atlantic the 
average storm track was displaced to the north in a similar 
manner for a similar reason. 

Storminess east of the Rockies reached its highest 
frequency during the last half of February. Essentially, 
there were two major storm tracks—one along the 
Canadian border, the other from the Southwest to the 
Northeast. Both tracks can be reasonably ‘‘fitted”’ to 
the flow at 700 mb. and to the mean sea level chart. 


Along the sea level track, associated with the new trough 
in the Plains, was found the most intense storm of the 
winter season in the East, as well as those systems re- 
sponsible for record-breaking precipitation in the South- 


east. 
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